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ABSTRACT 



The RAdial Velocity Experiment (RAVE) is a spectroscopic survey of the Milky Way which already collected over 400 000 spectra of 
~ 330 000 different stars. We use the subsample of spectra with spectroscopically determined values of stellar parameters to determine 
^^ , the distances to these stars. The list currently contains 235 064 high quality spectra which show no peculiarities and belong to 210 872 

^+ ■ different stars. The numbers will grow as the RAVE survey progresses. The public version of the catalog will be made available 

• ' through the CDS services along with the ongoing RAVE public data releases. 

r^ I The distances are determined with a method based on the work by Breddels et al. (2010). Here we assume that the star undergoes a 

^^ , standard stellar evolution and that its spectrum shows no peculiarities. The refinements include: the use of either of the three isochrone 

>— ^ ■ sets, a better account of the stellar ages and masses, use of more realistic errors of stellar parameter values, and application to a larger 

dataset. The derived distances of both dwarfs and giants match within ~ 21% to the astrometric distances of Hipparcos stars and 
to the distances of observed members of open and globular clusters. Multiple observations of a fraction of RAVE stars show that 
. ^ , repeatability of the derived distances is even better, with half of the objects showing a distance scatter of i 11%. 

ji^ ■ RAVE dwarfs are ~ 300 pc from the Sun, and giants are at distances of 1 to 2 kpc, and up to 10 kpc. This places the RAVE dataset 

}_( ' between the more local Geneva-Copenhagen survey and the more distant and fainter SDSS sample. As such it is ideal to address some 

C^ I of the fundamental questions of Galactic structure and evolution in the pre-Gaia era. Individual applications are left to separate papers, 

here we show that the full 6-dimensional information on position and velocity is accurate enough to discuss the vertical structure and 
kinematic properties of the thin and thick disks. 

Key words. Stars: distances - Catalogs - Surveys - Galaxy: structure - Galaxy: stellar content - Stars: statistics 

1. Introduction nates and in our case proper motions. But characterizing the 

three-dimensional position requires knowledge of the distance 

It is now widely accepted that our Galaxy can be used to the star and describing its motion requires knowledge of its 

as a testbed to study galaxy formation and evolution radial velocity. So one needs to measure also these quantities in 

dFreeman & Bland-Hawthorn] |2002|). A vital tool for this goal order to have complete information of an object's place in the 

is the ability to use complete information on position, motion 6-dimensional phase space which combines the position and ve- 

and physical properties of a large number of individual stars, lodty vectors. RAVE (RAdial Velocity Experiment) is the only 
Some of the quantities are usually known, i.e. the sky coordi- 
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large scale survey measuring radial velocities spectroscopically, 
with the exception of the SDSS survey which mainly samples the 
halo of the Galaxy. RAVE already collected over 400 000 spectra 
of ~ 330000 different stars. Most of the RAVE objects are thin 
and thick disk stars. So the SDSS and RAVE surveys are comple- 
mentary in the sample surveyed (halo vs. disk), magnitude range 
(fainter vs. brighter) and even the celestial hemisphere (North vs. 
South). The RAVE survey is described i n the articles accompa- 
nying the first two public data releases dSteinmetz et al. Il2006t 
IZ witter et al. 1120081 hereafter Z08). The article de scribing the 
third public release is in preparation (Sieb ert et al. 1 12010). With 
the exception of the first data release the subsequent ones also 
include results of spectroscopic measurements of stellar phys- 
ical parameters: the efifective temperature, surface gravity, and 
metallicity. 

The RA VE dataset can b e used to study a range of scien- 
tific topics dSteinmetz I l2003h . including present structure and 
formation of the Galaxy, remnants of recent accretion events, 
and discovery of individual peculiar objects and spectroscopic 
binary stars (Matiievic et alj l2010h . It was already used to 
better constrain the Galactic escape sp eed at the Solar radius 
(ISmith et alTI l2007h . IVeltz et al. I (l2008l) studied kinematics to- 
wards the Galactic poles and identified disconti nuities that sepa- 
rate t hin disk, thick disk and a hotter component. lSeabroke et al. I 
(12008 ) searched for in-falling stellar streams on to the local 
Milky Way, and ISiebert et aU ([2008) estimated the tilt of the 
stellar velocity ellip soid and its implications for mass models. 
Pasetto et al. ' ('2010') studied the kinematics of the thick disk. 
Munari et al. (2008) showed that a diffuse interstellar band at 
8620 A can be used to reliably measure the amount of galactic 
interstellar reddening. Munari et al. (2009) studied the proper- 
ties of Luminous Blue Variables in the Large Magellanic Cloud. 
They used a photo-ionization analysis of their rich emission line 
spectra to demonstrate the great diagnostic potential of RAVE in 
spectroscopy of peculiar stars. Additional ongoing studies have 
been listed in lSteinmetz et al.l (12006) . 

Most of these studies would profit from the knowledge of 
stellar distance. With very few exceptions, RAVE stars are too 
far away and too faint to have trigonometric parallaxes measured 
by the Hipparcos mission. The RAVE sample is magnitude lim- 
ited (9 < /denis < 13) but otherwise entirely randomly selected 
from the Tycho (bright stars), SuperCOSMOS (faint stars) and 
DENIS catalogs. So only a small fraction of objects observed 
by RAVE are members of stellar clusters or are of special type 
for which standard candle techniques can be used to judge their 
distances. One choice is to assume that all RAVE stars belong 
to the main seque nce and to use their photome tric colors to in- 
fer their distance (Kle ment. Fuch s & Rix 2008). Unfortunately, 
giants form about half of the RAVE sample, so any conclusions 
based on the assumption that they are main sequence stars will 
be biased. Another possibility is to use photometric colors to 
isolate the red clump stars. This is a viable technique, but it ap- 
plies only to a fraction of the sample and may be influenced by 
contaminating outliers of similar colors. 

Over 90% of the RAVE spectra belong to normal stars and 
have no substantial instrumental problems. So one may use the 
knowledge of spectroscopically measured physical parameters 
to infer their abolute position on the H-R diagram and thus 
their absolute magnitude. The choice of a near infrared passband 
(72MASs) and the fact that the vast majority of RAVE stars lie at 
high Galactic latitudes {\b\ > 20°) means that the influence of in- 
terstellar reddening is negligible. This approach was intr oduced 
by our first paper on the subject dBreddels et al. Il2010l here- 
after BIO). Here we build on this solid foundation and refine this 



approach by a better accounting for measurement errors, stellar 
evolution knowledge, allowing for three different sets of stellar 
isochrone computations, a refined computation method, and by 
expanding the sample to the recently observed stars. Another 
distance computation method using a Bayesian approach tested 
on a Galactic pseudo-dat a sample will also be use d to calculate 
distances to RAVE stars dBumett & Binnev"2010'). 

All these efforts are aimed at a common goal to derive the 
most reliable distances to RAVE stars. The goal is important 
and it is likely that different distance determination methods will 
be preferred for different types of stars. The situation could in- 
deed prove to resemble the problem of computation of stellar 
isochrones where it is best for the user to choose a particular set 
depending on the problem which is being considered. 

By the end of the decade a public catalog of extremely accu- 
rate and numerous astrometric distances measured by the ESA's 
mission Gaia is expected to become available. Still, as with any 
magnitude limited survey, there will be a significant portion of 
stars for which their distances are too large for accurate trigono- 
metric parallaxes to be determined. So photometric distance de- 
termination methods will remain to be useful also in the post- 
Gaia era of Galactic astronomy. 

The structure of the paper is as follows: Sec. |2] introduces 
the refined distance determination method. In Sec.[3]we test the 
results by comparing them to the available Hipparcos trigono- 
metric parallax measurements and to the distances of open clus- 
ters with some confirmed members observed by RAVE. Sec. |4] 
discusses statistical properties of the results and the structure of 
the publicly available catalog, to be followed by general conclu- 
sions in Sec. |5] The reader is encouraged to consult BIO for a 
more complete discussion of different methods to measure stel- 
lar distances and some of the scientific results. Detailed scien- 
tific results related to the structure and evolution of the Galaxy 
as well as distances of individual stars and stellar clusters will be 
pubUshed separately. 

2. Refined distance determination method 

BIO used val ues of stellar parameters and the Yonsei-Yale 
dDemarque et al. 2004) isochrones to derive the absolute /aiviASS 
magnitudes. Comparison with the apparent 2MASS J magni- 
tudes then yielded distances. The values of stellar parameters 
measured by RAVE are not error- free. BIO assumed constant er- 
rors in temperature (300 K), gravity (0.3 dex), and metallicity 
(0.25 dex) for all stars. So the values derived by RAVE were jit- 
tered according to these assumed standard uncertainties. 5000 
realisations of each observation were calculated by sampling 
Gaussian distributions in each observable. The final absolute J 
magnitude was the mean of the absolute J magnitudes pertaining 
to the closest isochrone matches for each of the realisations. The 
closest match was derived by minimizing the usual ;^'^ statistics. 
The distribution of absolute J magnitudes for individual realisa- 
tions also gave the standard error of the derived absolute magni- 
tude and thus of the distance to the star. 

In formulating refinements to the BIO approach we alter 
the method for distance derivation. Instead of using the closest 
match, an expectation value of the absolute magnitude is derived 
considering several isochrones at a given metallicity. Our goal is 
to derive the most likely distance to the star So we assume that 
we are dealing with a normal star that follows stellar evolution 
as mirrored by theoretical isochrones. We use our knowledge of 
stellar evolution to compensate for some of the errors in values 
of stellar parameters as derived by RAVE. So the results of this 
paper should be very useful in Galactic and general population 
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studies. But we note that our results may not be applicable to 
peculiar objects or to those in rapid transit phases of stellar evo- 
lution. The latter because the adjacent longer lived and so more 
populated evolutionary phases may be within the measurement 
errors. The refinements to the BIO method are described below. 



2. 1 . The three sets of isochrones 

The distance c omputation of BIO was based on the Yonsei-Yale 
(hereafter YY, iDemarque et al. I '2004) set of isochrones. Here 
we add two additional sets: the Dartmouth (iDotter et al. Il2008h 
and the Padova isochrones jBertelli et al. Il2008h . While the for- 
mer two sets are very useful for the derivation of properties of 
stars which are on or close to the main sequence, their com- 
putations do not extend to the asymptotic giant branch (AGB) 
or red clump (RC) evolutionary phases. The Padova set on the 
other hand nicely covers both the AGB and the RC phase. So it 
is particularly useful to derive distances to giant stars. All sets 
of isochrones tabulate also the absolute J magnitude. In the case 
of YY and Dartmouth sets the absolute J magnitude can be di- 
rectly compared to the a pparent one measu red by the 2MASS 
survey, while eq. 1 from iKoen et al. I (l2007h was used to trans- 
form the Padova set into the 2MASS system. In all cases we 
used isochrones without any enhancement of alpha elements. 
This quantity is not measurable reliably by RAVE and has a 
very small influence on the final distance results (see BIO for 
a detailed discussion). 

All isochrone sets cover the parameter space in the neces- 
sary detail and in some cases even ofi'er interpolation routines 
to derive values at points in between. These routines were used 
when possible, otherwise a standard linear interpolation scheme 
was used. So we generated a fine grid of isochrones with age, 
mass and metallicity steps much smaller than the observational 
errors. The mass step was always smaller than 0.01 M©, and the 
metallicity step was 0.1 dex. Sampling in age is discussed below. 

All computations were performed three times, deriving a 
separate set of distances for each set of isochrones. We publish 
all three distances and their errors. It is then up to the user to 
decide which one is the most suitable for the type of problem or 
the type of object that is being considered. Details of isochrone 
computation and the corresponding assumptions are described 
in detail in the original papers presenting each isochrone set. 



2.2. Isochrones uniformly spaced in time 

BIO used a set of isochrones with 40 diff'erent ages, spaced loga- 
rithmically between 0.01 and 15.0 Gyr. This implies a time step 
of only 2 Myr at the young end and up to 2.5 Gyr at the old 
end of the isochrone set. Such a choice covers the [temperature, 
gravity] plane well. But the use of the closest match implies an 
implausibly high probability that a certain realisation lies within 
a rapid evolutionary phase, such as a contraction towards the 
main sequence or a transit from the main sequence towards the 
region of giants. The bulk of the Galactic populations is expected 
to be found in slow, long-lasting phases of stellar evolution. So, 
although it may be possible that RAVE observed some objects 
in rapid phases of stellar evolution, this is not probable in indi- 
vidual cases. 

In the present calculations the isochrones are spaced uni- 
formly and not logarithmically in time. So the density of grid 
points is inversely proportional to the speed of motion of a star 
across the [temperature, gravity] plane. We use a step of 20 Myr 
for ages up to 12 Gyr. The youngest age depends on its avail- 



ability for a particular isochrone set. It is 0. 1 Gyr for the YY set, 
1 Gyr for Dartmouth and 0.5 Gyr for the Padova one. Uniform 
spacing in time gives a proper weight to longer evolutionary 
phases, while rapid transit periods have their importance dimin- 
ished. The adopted age range includes the span of ages of a gen- 
eral stell ar population as determined by the Genev a-Copenhagen 
survey (iHolmberg. Nordstrom & Andersen II2009I) . 

2.3. Weighting according to mass and absolute magnitude 

The quantity changing along a theoretical isochrone is stellar 
mass. Isochrones are tabulated so that connecting adjacent points 
follows their shape well. So the points are tabulated in relatively 
large mass steps while the isochrone lies on the main sequence. 
The mass steps get much smaller when the isochrone bends to- 
ward giants. 

We use a scheme that weights individual solutions. Clearly 
the weight of a given point should be proportional to a mass 
step dm which is an average of mass differences to the adjacent 
mass points on a given isochrone. It should also be noted that 
high mass stars are much less frequent than ones with sub-solar 
mas ses. dN/dm is the Galactic mass function as given in eq. 17 
of Chab rier I (l2003h for sub-solar stars and in his eq. 1 1 for stars 
more massive than the Sun. Moreover, more luminous stars at a 
given apparent magnitude are farther away than their less lumi- 
nous counterparts. The volume of space within which the stars 
with a given absolute magnitude M/ in the Cousin / band appear 
brighter than the given apparent limiting magnitude is propor- 
tional to lO""'^*''. We use the / band, as both the surveys used to 
select the observing sample and the final RAVE spectra are cen- 
tered on this band. If we assume that the space density of stars 
does not depend on distance the following relative weight x of 
points along the isochrone is obtained: 



X — dm(dN/dm)lO 



-0.6M, 



(1) 



One may argue that the dN/dm term can be used only for 
a volume limited survey, while RAVE is a magnitude limited 
one. Similarly, the stellar distribution is not uniform. Stars at 
large distances are quite far away from the Galactic plane, so the 
frequency of these luminous, massive stars is even lower than 
would be in the case of a uniform space density. Fortunately 
these assumptions are not crucial. It is reasonable to take into 
account that massive stars are uncommon and that RAVE is able 
to sample a larger volume of more luminous stars. But the ex- 
act mathematical form of these prescriptions is not so important, 
and we can choose the simplest relation using varying weights. 
The reason lies in the quality of RAVE's measurement of values 
of stellar parameters. Their errors are rather small, even though 
not negligible. So any large offsets in temperature, gravity and 
metallicity are virtually excluded. The role of equation[T]is to as- 
cribe a relative weight to points on the isochrones that are close 
to the initial RAVE values. 

As a test we recalculated the whole database using x = dm. 
So we ignored any mass or luminosity dependence. The results 
are illustrated in Fig. [1] We see that most stars are closer than 
in the case when we use eq. [T] The magnitude term that favours 
more luminous and hence more distant stars in eq. [1] therefore 
dominates over the mass function dN/dm which disfavours mas- 
sive (and so luminous) stars. But differences in distances derived 
by the test procedure and the usual one are small: the mode of 
the distance ratio distribution derived by the two methods is 0.94, 
and the median is 0.87. So even the extreme simplification which 
ignores any assumptions used to derive eq. [T| changes the dis- 
tances by only ~ 10%. 
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Fig. 1. Ratio of distances derived under different assumptions. 
di is derived assuming that the probability of a star to be ob- 
served does not depend on its mass or brightness, while the dis- 
tance d2 takes into account the mass and brightness functions, as 
described in the text. 



2.4. Assumed errors of stellar parameters 

The errors of stellar parameter values measured by RAVE de- 
pend on stellar type. Typical temperature errors for stars cooler 
than 9000 K are around 400 K. The errors are smallest for su- 
pergiants. Their atmospheres are the most transparent ones, so 
the wealth of spectral lines arising at different optical depths can 
improve the temperature accuracy. Understandably the errors for 
metal-poor stars are larger than for their Solar counterparts. The 
errors get considerably worse for stars hotter than 9000 K, where 
most metal lines are missing and the spectrum is largely domi- 
nated by hydrogen lines. Similarly, strong wings of hydrogen 
lines, which are sensitive to gravity, allow small gravity errors in 
hot stars. On the other hand rather narrow metallic lines in the 
RAVE wavelength range, including those of Ca II, do not allow 
an accurate determination of gravity in cool stars. The gravity 
error in cool stars has a strong gravity dependence: in dwarfs 
it is ~ 0.6 dex, but the rather transparent atmospheres of giant 
stars still allow for a more accurate gravity determination, with 
a scatter that can be as low as 0.2 dex. Finally, the typical error 
in metallicity for stars cooler than 7000 K is 0.2 dex. The error 
for hotter stars is understandably much larger, as these stars lack 
most of the metallic lines in their spectra. 

Given this vivid behaviour it seems important to allow for a 
variation of the accuracy of stellar parameter values across the 
parameter space. This is a refinement over other studies aim- 
ing at distance computation using values of stellar parameters as 
determined by RAVE. Here the errors are obtained from a lin- 
ear interpolation of results presented in Fig. 19 and equations 
22 and 23 of Z08. Note that these er ror estimates may be rather 
conservative (ISteinmetz et al. Il2009h . 

2.5. Calculating the distance modulus 

Values of stellar parameters, as determined by RAVE, are not 
error- free. So BIO first determine the closest match to the RAVE 
values in the isochrone parameter space. This closest match is 
then jittered according to the assumed parameter uncertainties. 



Each of the realizations again yields the closest match in the 
isochrone space, and a suitable average of these closest matches 
yields the final value of the distance and its scatter 

This procedure can benefit from a refinement regarding the 
use of the closest match. Imagine we have a Solar type main- 
sequence star that would have the reported values of RAVE pa- 
rameters placing it just above the main sequence of the H-R dia- 
gram. Jittering the reported value of its gravity would mean that 
half of the realizations would imply a higher gravity and so a 
main-sequence solution. But the other half of the jittered real- 
izations would claim the star is above the main sequence in a 
rapid transition towards the giant phase. The use of the closest 
matches would thus produce half of the realizations involving 
massive stars in a rapid transit toward giants and half of main- 
sequence cases. The relative weight of realizations involving the 
rapid transition phase is therefore too large. It also influences 
the value of the mean of absolute / magnitudes for individual 
realizations and so the derived distance to the star 

We use a modified procedure to calculate the most likely 
value of the absolute J magnitude. The distance modulus then 
follows from its difference with the apparent J magnitude as 
measured by the 2MASS survey. We do not jitter the mea- 
sured parameter values and calculate individual realizations to 
be solved with the closest match. Instead, the expectation value 
of the absolute magnitude is derived using the relative weights 



X(a, m) — x{a, m) e 



-y((3,m) 



(2) 



[r(RAVE)-r(fl,m)]^ [log^(RAVE)-logg(fl,m)]2 

y(a, m) — ■ — H 



^" r(RAVE) 



^"logg(RAVE) 



where a is age and m mass at a given point on the isochrone 
with the metallicity value of the RAVE measurement. Standard 
deviations of the measured value of temperature io-T(RAVE)) and 
gravity (criog^(RAVE)) are determined as described in Sec. 12. 41 The 
final value of the absolute J magnitude (Mj) and its standard 
deviation {cr(Mj)) are calculated from the first moments of the 
distribution 



Mj 



I^a Zm Mj(a, m)X(a, m) 



cr'(Mj) 



Za Zm X(a, m) 

ZaZ,nX(a,m)[Mj - Mj{a,m)f 
Za Zm X{a, m) 



(3) 



(4) 



The procedure makes a brute force average of all solutions for 
an isochrone set of a metallicity measured by RAVE. Typical 
values of metallicity errors are ^ 0.2 dex, or about half of the 
step in a grid of stellar atmosphere models used by Z08 (their 
Fig. 19). On the other hand the errors in temperature or grav- 
ity are typically larger than the grid step. The isochrone sets for 
different values of metallicity are very similar, as already dis- 
cussed by BIO (see its Fig. 2). A lower metallicity corresponds 
to a more transparent stellar atmosphere with a hotter eff'ective 
temperature, but has little effect on the stellar energy output. So a 
metallicity change causes a shift of the isochrone in the tempera- 
ture direction, but the absolute magnitude stays very similar (see 
isochrone comparison in Fig. |2j5 below). It is thus not surpriz- 
ing that the results on absolute magnitude show no significant 
changes when using metallicity as another independent variable. 
In particular, the distances do not change by more than 10% for 
90% of the spectra when using metallicity as an independent 
variable. But an extra term significantly prolongs the calcula- 
tion. So we decide to adopt the simplified form given in Eq. |2] 
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If isochrones with the value of RAVE-measured metallicity do 
not exist they are obtained by a linear interpolation from the ad- 
jacent ones. Similarly, some of the isochrone sets are calculated 
for a set of iron abundances ([Fe/H]) which in general differ from 
the metallicities ([M/H]) measured by RAVE. Equation 21 from 
Z08 is used to convert between iron abundance and metallicity. 
Similarly we use eq. 20 from Z08 to calibrate the metallicity 
values. 

We explained in Sec. 12.21 that the isochrones are uniformly 
spaced in age, while Sec. |23] discusses relative weights x(a,m) 
along an isochrone. A Gaussian form of eq. |2] means that only 
points on the isochrones that have the temperature and grav- 
ity within a few standard deviations from the RAVE measure- 
ment are important. This scheme weights each solution accord- 
ing to its distance from the position measured by RAVE. So the 
weighting goes by distance in the [temperature, gravity] parame- 
ter space and not by area, as is the case with the closest match. So 
we avoid giving excessive weights to isochrones in the scarcely 
populated regions of the parameter space, where a given solution 
is the closest match to a substantial range of the [temperature, 
gravity] pairs. 



3. Tests of the derived distances 

Stellar distances are generally not known, so we can test our 
results only for certain subsets of the measured sample. Two 
such cases will be considered: Hipparcos stars and confirmed 
members of stellar clusters. In future we plan to test the result s 
also with a Galactic pseudo-data set jBurnett & Binnev1l2010h . 
Combined results of these tests and experience obtained through 
a scientific use of distance results will be discussed in future re- 
leases of RAVE data. 



3.1. Hipparcos stars 

Hipparcos stars h ave their distances m easured through trigono- 
metric parallaxes jvan Leeuwen 112007 ). Our goal is to use ob- 
jects with well determined trigonometric distances and with 
high-quality RAVE spectra. This choice can then allow us to 
check for systematic errors in our distance computation which 
is the primary goal of our test. The formal errors on the paral- 
laxes are 1 .0 + 0.5 mas, so we adopted a lower limit of 7.7 mas 
for the parallax, corresponding to a distance limit of 130 pc. 
Moreover, for this test we considered only RAVE spectra with 
a S/N ratio > 40. RAVE observed 1525 such stars so far. Finally 
we require that a formal error in spectroscopic distance modulus 
< 0.6 mag (corresponding to a distance error of 32%). The last 
constraint is satisfied by 725 spectra for calculations based on 
the YY isochrones, 776 spectra for the Dartmouth set, and 674 
for the Padova set. 277 of them had their spectroscopic distances 
already included in BIO. 

RAVE confirmed that most of these stars have a Solar-like 
composition ([M/H]= -0.12 + 0.25). Fig. |2}i shows RAVE mea- 
surements of the effective temperature and surface gravity. As 
expected nearly all stars are dwarfs. The magnitude limits of the 
RAVE survey also imply that most stars are somewhat hotter 
(and brighter) than the Sun, lying in the general turn-off region. 
As discussed in Sec. l2.4l their temperature determination is rather 
accurate, but their gravity errors are quite high. So RAVE is able 
to determine their general main-sequence type, but not any fur- 
ther details. The error on gravity is even worse for cool stars, 
so the small group of deviant stars at the bottom-right corner of 
Fig. |2^ can easily be normal main-sequence stars. 
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Fig. 2. HR diagram of Hipparcos stars closer than 130 pc. The 
top panel (a) shows RAVE measurements of gravity and effec- 
tive temperature, while the bottom panel (b) presents the result 
of our calculation of the absolute J magnitude using the YY set 
of isochrones. Grey points are Hipparcos stars with a trigono- 
metric distance error smaller than 20%, a spectroscopic distance 
error below 32% and with RAVE spectra with S/N > 40. Black 
points are the subset where the BIO and the three isochrone set 
calculations presented in this paper found a spectroscopic dis- 
tance error smaller than 32%. Black and gray lines trace YY 
isochrones with Solar and subsolar composition and for ages be- 
tween 1 and 1 1 Gyr with a 1 Gyr spacing. The vast majority 
of Hipparcos stars close enough to have accurate trigonometric 
parallaxes are dwarfs. 

The top panel of Fig. |2] presents the input for the spectro- 
scopic distance calculation, while the bottom panel depicts its 
output. The ordinate is now the calculated absolute J magni- 
tude resulting from the YY set of isochrones. We see that the 
method nicely corrects the too-high gravities from the top panel 
and places the sub-dwarf stars back on the main sequence. Also 
the cool stars supposedly lying above the main sequence are now 
nicely settled onto it. 

Fig. |3] compares trigonometric distances derived by 
Hipparcos (dj^) to the RAVE spectroscopic distances (d) as de- 
termined by BIO and by the present calculation using the three 
isochrone sets. Our restriction to very high-quality measure- 



Zwitter et al.: Distance determination for RAVE stars using stellar models, n 



Table 1. Comparison of spectroscopic and trigonometric distances 



(1) 



Calculation 



Sample A*'* 
k a spectra 



Sample B'*' 
k a spectra 



BIO 

This paper (YY isochrones) 

This paper (Dartmouth isochrones) 

This paper (Padova isochrones) 

Notes. (1) Spectroscopic (d) and trigonometric {d^) distances are assummed to follow the relation d = kd^. (2) Hipparcos stars with a trigono- 
metric distance error smaller than 20%, a spectroscopic distance error below 32% and with RAVE spectra with S/N > 40. (3) A subset of Sample 
A where all four calculations find a spectroscopic distance error smaller than 32%. 



0.92 ±0.01 0.18 277 0.89 ±0.01 0.18 207 

1.03 ±0.01 0.23 725 1.09 ±0.02 0.25 207 

1.01 ±0.01 0.23 776 1.07 ±0.02 0.24 207 

1.03 ±0.01 0.25 674 1.09 ±0.02 0.26 207 




50 100 

d, [pel 

Fig. 3. Comparison of distances to Hipparcos stars. The trigono- 
metric parallax distances dj^ are compared to spectroscopic dis- 
tances d from BIO (a), and to values derived here using the YY 
(b), Dartmouth (c), and Padova (d) isochrones. We are using the 
same selection criteria and color scheme as in Fig.|2] The error- 
bars show average formal errors for stars in four distance bins. 
The diagonal line is a 1 : 1 relation. 





Fig. 4. Distances to open clusters (^/cluster) from Dias et al. 
(2002, 2007) compared to spectroscopic distances d of cluster 
members observed by RAVE from BIO (a), and to values de- 
rived here using the YY (b), Dartmouth (c), and Padova (d) 
isochrones. Circles denote dwarfs (loggRAVE > 3.8) and trian- 
gles likely giants (loggRAVE ^ 3.8). Open symbols mark the 
likely members of M 67. The line marks the 1:1 relation. The 
pairs of distances in the shaded area are equal to within 20%. 



ments of trigonometric parallax allows us to ignore any trans- 
formation bias due to the non-linear relation between parallax 
and distance. 

The spectroscopic distance determinations are generally in 
good agreement with their trigonometric parallax values. Table[T] 
reports the best fits of the relation d - kd^ for two samples and 
for four different sets of calculations. The slopes k are close to 
unity. The BIO distances are somewhat underestimated, while 
the ones reported in this paper are slightly overestimated if 
compared to trigonometric parallax measurements. On the other 
hand the scatter in fc is a bit larger for measurements in this pa- 
per than was the case in BIO. This is partly due to a difference 
in the assumed errors of stellar parameters. While BIO use the 
same error for all spectra, our errors are variable and depend on 
the type of star and the S/N ratio of the spectrum. Our assumed 
errors are generally larger, a typical error in gravity is 0.45 dex 
compared to 0.3 dex used by BIO. Assuming fixed error values 
from BIO (300 K in temperature and 0.3 dex in gravity) we get 
more satisfying fits for the Sample B: fe = 1.03 ± 0.01, cr = 0.21 
(YY isochrones), k = 1.04 ± 0.01, cr = 0.21 (Dartmouth), and 
k= 1.05 ± 0.01, cr = 0.22 (Padova set). 

Note that many of the Hipparcos stars with the temperature 
close to 6000 K have log g ~ 4.8, a value too high for this part of 



the main sequence (Fig.|2^). This offset challenges the distance 
routine which assumes a Gaussian error distribution. The offset 
is more frequent for the Hipparcos stars than for typical RAVE 
spectra. About 15% of the Hipparcos stars hotter than 5500 K 
have a claimed \ogg > 4.8, compared to just 8% of the cases 
for the general RAVE sample. Also, Hipparcos stars sample the 
cool part of the main sequence down to temperatures of 3500 K. 
But the general RAVE sample has hardly any dwarfs cooler than 
4500 K, with a vast majority at around 6000 K. The difference is 
due to a magnitude limited nature of the survey, where cool stars 
are sampled only in a small volume around the Sun. 

So the scatter obtained for the Hipparcos stars may not be 
representative of the entire population of RAVE spectra of the 
main-sequence stars. Note that the Hipparcos stars are generally 
brighter than the general sample observed by RAVE, so their 
spectra have a higher S/N ratio. Thus the distance error esti- 
mates for the Hipparcos stars may be somewhat optimistic. It 
is interesting to note that the smallest distance errors are fre- 
quently found for stars with overestimated gravity values. This 
is understandable, as for such stars the main-sequence solution 
is an obvious choice. On the other hand the stars with more real- 
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istic and so lower values of gravity could also be turn-off stars, 
so their absolute magnitude determination is less certain. 

To summarize, trigonometric parallax measurements and re- 
sults of our spectroscopic distance computations show average 
discrepancies of .^5%. The scatter for measurements of individ- 
ual stars is i 24%. Including the formal error of the Hipparcos 
parallaxes in quadrature, we find that the spectroscopic distances 
are accurate to within i 22%. 



3.2. Distances to members of star clusters 

RAVE observed virtually no giant stars with an accurate mea- 
surement of their trigonometric parallax. Our spectroscopic dis- 
tance measurement, however, can easily reach the giants that 
form about half of the RAVE sample. So we used known mem- 
bers of open and globular clusters observed by RAVE to extend 
our check of the distance scale. 

The Vizier catalog of open clusters of iDias et al. I (l2002h . 
Version 2.8 was used. A star is considered to be a known open 
cluster member if: (i) it lies within the apparent diameter of the 
cluster center in both coordinates , (ii) its pro per motion and the 
value for the cluster reported by iDias et alTl (12002) do not dif- 
fer by more than the sum of their reported uncertainties, (iii) 
the same is true also for radial velocities (if iDias et al. 1 12002 , 



do not report the radial velocity error a difference of 2 km s 
is considered acceptable), (iv) the Simbad catalog reports that 
the object is a "star in cluster" (class *iC). Clearly these re- 
quirements are conservative. They exclude many potential clus- 
ter members observed by RAVE. But our goal here is to use defi- 
nite clustermembership to check our distances. We also checked 
the [Harris et al. I ([1996) catalog of globular clusters. The star 
J 125905. 2-705454 is 3.43 arcmin from the center of the glob- 
ular cluster NGC 4833, which has a half-mass radius of 2.41 ar- 
cmin. Its heliocentric radial velocity of 202.1 ± 0.8 km s"' (as 
measured by RAVE) matches the cluster's radial velocity of 
200.2 ± 1 .2 km S-' . Similarly its RAVE metallicity [M/H]= -1 .2 
is consistent with the cluster iron abundance of [Fe/H]= -1.79. 
So we consider it to be a member of this globular cluster No 
other members of globular clusters were found. 

Fig.|4]and Table[2]summarize the results. Altogether we have 
observations of 11 secure members of 7 different open clus- 
ters, one globular cluster member and 3 additional likely cluster 
members. Four stars are dwarfs (logg > 3.8), the rest are gi- 
ants (iogg < 3.8). Our results show a satisfactory performance 
for dwarfs and for giants. Globular cluster member J125905.2- 
705454 is omitted from Fig. |4] due to scale. Table |2] shows that 
also for this supergiant (Teff - 4250 K, log g - 0.8) the spectro- 
scopic distance is consistent with the value from the literature. 
The derived spectroscopic distances are generally within 20% 
of the distance to th e cluster as given in iDias et alTI ( |2002|) and 
[Harris et al. I (Il996l) . The average and the scatter of the spectro- 
scopic to cluster distance ratio for the entries in Table |2] equals 
0.97+0.14 (1.11 ±0.08), 1.03+0.15 (1.16±0.08), and 0.89±0.15 
( 1 . 1 1 + 0.08) if using the Y Y, Dartmouth, and Padova isochrones 
respectively. The values in brackets are obtained if we fit only 
the open cluster members. The difference of solutions inside and 
outside the brackets is quite large, given that it corresponds to a 
contribution of a single globular cluster member So we note that 
all values are consistent with the 1:1 relation to within ~ 10%. 

BIO use observations of the open cluster NGC 2682 (M 67) 
to study the distance scale of giants. The distances found are 
generally too large. Here we show that our distances are consis- 
tent with the adopted value of 910 pc from the literature. The 
two members in Table |2] give 930 ± 160 pc (YY isochrones). 



980 + 160 pc (Dartmouth), and 880 + 150 pc (Padova). If we 
would add the likely M 67 members by their heliocentric ra- 
dial velocity (31.2 km s"'< RV < 33.4 km s"') the number of 
observed M 67 members would increase to five and the results 
would be 930 ± 110 pc (YY), 980 + 110 pc (Dartmouth), and 
940 ± 100 pc (Padova). 

Stellar clusters have a known value of the color excess Eb-v, 
so they allow some discussion of the influence of interste l lar ex- 
tinction on our distance estimates. iFiorucci & Munaril (l2003h 
give for the extinction in the J band Aj ^ 0.887£'b_v'. Nearly 
all clusters in Table|2]have Eb-v < 0.1, so their extinction in the 
J band is Aj < 0.089. Ignorance of extinction makes their dis- 
tances overestimated by no more than 4%. The only clusters with 
the color excess larger than 0. 1 are NGC 2477 ( Eb-v = 0.24, 
'Dia setal.ll2002l) and NGC 4833 (Eb-v = 0.33, ' Harris et al. I 
19961). Neglecting interstellar extinction should thus make their 
distances overestimated by 10% and 14%, respectively. Note that 
these two clusters lie very close to the Galactic plane (b - -5.8° 
and b = -8.0°, respectively). A vast majority of the RAVE sam- 
ple has \b\ > 20°, so its interstellar extinction can be safely ig- 
nored. 

RAVE observed more objects that are likely members of 
open clusters and satisfy the same selection criteria, with the 
difference that Simbad does not recognize them as cluster mem- 
bers. There are a number of such likely members of NGC 2423, 
and a few objects in the open clusters CoUinder 121, ASCC 33, 
and Alessi 5. Here we compared spectroscopic distance mea- 
surement of confirmed cluster members to known cluster dis- 
tances. A study of possible cluster members and applications us- 
ing spectroscopic distances to stellar clusters will be published 
elsewhere. 



4. Derived distances 

The method described in Sec. |2] has been applied to a sample 
of RAVE stars with spectroscopically determined values of tem- 
perature, gravity and metallicity. First we describe the sample 
and stellar properties as reflected in the derived absolute magni- 
tudes. Then we move to statistical properties of derived distance 
distributions, including their error estimates. We conclude with 
a detailed discussion of the catalog and its availability. 



4. 1 . The sample 

RAVE is an ongoing survey that uses the multi-fiber spectro- 
graph of the UK Schmidt telescope at the Anglo Australian 
Observatory to obtain stellar spectra with a resolving power of 
~ 7500 in the wavelength range from 8410 A to 8795 A. It is 
conducted by an international collaboration led by M. Steinmetz 
and aims to obtain a million stellar spectra in the magnitude 
range 9 < /denis < 13. The mode of the signal-to-noise ratio 
per pixel of the > 400 000 that have so far been gathered follows 
the relation S/N = io-0-2(^denis-i9.1) (^08). A detailed description 
of the survey is giv en in the papers accompanying the first tw o 
public data releases dSteinmetz et al. 12 006*: Zwitter et al. l 2008l) . 
soon to be followed by the third one ([Siebert et al. 2010). 

This paper is based on the sample of 332 747 currently re- 
duced spectra with a measurement of radial velocity. For 292 261 
of these we have determined the values of stellar parameters, 
with the rest being either too noisy or suffering from systematic 
problems such as a contamination by the second interference or- 
der A subsample of 25 1 364 have S/N > 20, so dependable es- 
timates of the values of stellar parameters (temperature, gravity 
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Table 2. Distances to members of star clusters. 



Object-ID 


RAj2000.0 


Decj2ooo.o 


Cluster 




S/N 


Spectroscopic Distances from RAVE 


![pc] 




[deg] 


[deg] 


ID 


d[pc] 




BIO 


Yonsei-Yale 


Dartmouth 


Padova 


J000128.6-301221 


0.369417 


-30.205861 


Blanco 1 


270 


64 


260±80 


270+70 


260+60 


260+60 


J000324.3-294849 


0.851583 


-29.813722 


Blanco 1 


270 


64 


200±30 


280+50 


270+50 


280+50 


OCL00148.1373319 


114.302500 


-13.748667 


NGC 2423 


770 


73 


2150±960 


960+220 


1040+220 


1100+190 


OCL00147_1373471 


114.434792 


-14.837000 


NGC 2422 


490 


44 


890±400 


610+90 


620+90 


710+160 


J075214.8-383848 


118.061750 


-38.646889 


NGC 2477 


1300 


50 


2110±1200 


1280+310 


1330+310 


1160+320 


J075242.7-382906 


118.178125 


-38.485083 


NGC 2477 


1300 


61 


2600±1300 


1690+320 


1770+320 


1720+340 


M67-0105 


132.821250 


11.804444 


NGC 2682 


910 


89 




800+170 


840+170 


740+150 


M67-0135(*) 


132.839833 


11.768361 


NGC 2682 


910 


50 




930+280 


1000+290 


1120+280 


M67-0223 


132.932875 


11.945083 


NGC 2682 


910 


57 


2510+700 


1070+260 


1130+270 


1010+250 


M67-2152(*) 


133.045750 


11.530361 


NGC 2682 


910 


48 


1900±940 


1000+270 


1040+270 


960+200 


M67-6515(*) 


133.069083 


11.327361 


NGC 2682 


910 


99 


1540±770 


860+170 


910+170 


880+190 


OCL00277_2236411 


165.825000 


_58.488444 


NGC 3532 


490 


52 


640±270 


610+120 


590+120 


580+120 


OCL00277_2236511 


165.934708 


-58.960972 


NGC 3532 


490 


39 


560+250 


470+150 


440+130 


450+150 


T775 1.00502.1 


171.452250 


-43.164611 


NGC 3680 


938 


67 


2670+990 


1210+280 


1270+280 


1140+260 


J 125905. 2-705454 


194.771667 


-70.914972 


NGC 4833 


5900 


37 




5500+1000 


5900+1100 


4900+1100 



Notes. (*) A likely cluster member: Simbad does not list it as a cluster member, but it is within the cluster area and has a radial velocity consistent 
with the cluster membership. 



and metallicity), which are virtually free from systematic off- 
sets, could be obtained (Z08). The S/N ratio depends only on 
the /denis magnitude of the object, on the observing conditions, 
and on the performance of the instrument. Objects to be observed 
are magnitude selected and observed in a random order. So our 
decision to discuss only objects with S/N > 20 is equivalent to a 
modified magnitude cut and can be taken into account when dis- 
cussing statistical properties of the sample. A few percent of the 
stars turn out to be double-lined spectroscopic binaries, coronal 
emission objects, or show other peculiarities. Our distance esti- 
mates assume the spectrum belongs to a single normal star. So 
these peculiar objects need to be e xcluded, for which we used 
an automated classification routine (^Matii evic et al. 1120101) . The 
routine unavoidably fails to identify some of the binary stars. In 
the case of a single-lined binary star the distance would be un- 
affected. On the other hand the distances to unresolved double- 
lined spectroscopic binaries are underestimated by up to a factor 
of \1. But such worst cases are rare (Ma tiievic et al. Il2010h . If 
a star in a binary is two magnitudes brighter than its unidentified 
companion, the error in its estimated distance drops to 8%. The 
final sample of stars for which the distances have been estimated 
contains 235 064 spectra with S/N > 20 and without detected 
pecuUarities. They belong to 210 872 different stars. 

4.2. Derived absolute magnitudes 

The immediate result of our method is an estimate of the ab- 
solute magnitude of the star. We use the 72mass passband be- 
cause the apparent magnitudes measured by the 2MASS survey 
are accurate, available for virtually all stars in our sample (Z08, 
Table 7) and sufficiently red to allow us to ignore any effects of 
interstellar extinction (see Sec. 13.2b . 

The properties of the sample can be studied with the H-R 
diagram in Fig.|5] which shows the results of the absolute mag- 
nitude calculations of BIO and by the method presented here for 
all three isochrone sets. The relative density of stars is repre- 
sented by linearly spaced colored contours. A subset of Solar 
metallicity isochrones is overlaid for illustration purposes. 

Two groups of stars are apparent: the moderate temperature 
stars which are close to the main sequence and the group of cool 
giants. Cool dwarfs are almost missing. This can be understood, 
because the RAVE survey has a fairly limited magnitude span 



(9 < /denis < 13), so the volume of space with cool dwarfs 
of appropriate apparent magnitude is smaller than the one with 
hotter dwarfs and much smaller than the volume occupied by lu- 
minous and thus distant giants. Hot dwarfs are also infrequent, 
which is not surprising given the low probability of their forma- 
tion and their short lives. Similar selection effects diminish also 
the relative frequency of supergiants. As a result of these con- 
siderations RAVE observes mostly dwarfs slightly hotter than 
the Sun and giants with surface temperatures of ~ 4700 K. 

Figure |5^ leaves some room for improvement. Many of 
the main-sequence stars seem to be just leaving the main se- 
quence, which is unlikely. Note that this is even more true if 
their metallicities were lower than Solar. Also, giants lie above 
the Solar-metallicity isochrones. This would be fine if they were 
very metal-poor, but RAVE metallicities tell that they are not. 
Consequently BIO find that the distance to M67 giants is 1480 
pc (if excluding the red clump stars), a value 60% larger than in 
the literature. The BIO distances to Hipparcos dwarfs however 
turn o ut to be in agreement w ith astrometric distance determina- 
tion of Ivan Leeuwen I (l2007h . 

Results using the refined method are given in Fig. |5j)-d. 
These panels show the results for the sample of 145 028 nor- 
mal spectra with S/N > 20 for which BIO also give distances 
(Fig.|5^). The whole sample we present in our Catalog (Sec. 14.6b 
is larger and includes 235 064 spectra. Both samples are drawn 
from the same selection strategy and observation protocol, im- 
plying very similar errors in derived values of stellar parame- 
ters and thus distances. Results of the refined computations show 
that dwarf stars now lie very close to the main sequence and the 
giants are in the area of moderately metal-poor giants with an 
age of a few Gyr. A closer inspection reveals that most main- 
sequence stars are more luminous than the locus of the zero-age 
main sequence. This is understandable, as a main-sequence star 
with ~ 6000 K and an age of ~ 5 Gyr is close to the turn-off point 
(see the 5 Gyr isochrone on Fig.|5]l, so it is notably brighter than 
it was when younger. 

Results for the calculations using the YY (Fig. |5j5) and 
Dartmouth (Fig. |5};) isochrones are very similar. However, the 
graph based on Padova isochrones (Fig.|5}l) shows a notable con- 
centration of giants at Mj ^ -0.9. This difference arises because 
only Padova isochrones include a proper treatment of the red- 
clump region. Since stars spend a long time in this region we 
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Fig. 5. Frequency of occurence of solutions in the absolute J magnitude vs. temperature diagram for four sets of calculations: (a) 
results from the BIO paper, (b-d) new calculations using the YY, Dartmouth and Padova isochrones, respectively. All panels include 
the same set of 145 028 spectra common to all four calculations, but the appearance of (b-d) for the whole dataset of 235 064 stars 
would be almost identical. Fifty colored iso-frequency contours are linearly spaced between 20 and 1900 stars per bin with a step of 
40. Each bin spans 50 K in the temperature and 0.2 mag in the magnitude direction. Gray lines trace Solar composition isochrones 
between 1 and 1 1 Gyr with a 1 Gyr spacing. 



expect many of the giants observed by RAVE to lie in the red 
clump. Hence, distances based on the Padova isochrones are to 
be preferred over ones calculated using the YY or Dartmouth 
isochrones, at least for giants. The locus of the red-clump stars 
lies at a bit hotter temperatures than the RGB stars with the same 
absolute magnitude. So the red-clump stars in Fig.|5}i are not su- 
perimposed on the red-giant branch. 

The red clump region in Fig. |5}l can be studied further. 
Fig. |6^ shows a histogram of Mj for temperatures between 
4500 K and 5000 K. Note the pronounced red clump maximum, 
modelled as a Gaussian with the peak at Mj - -0.91 ± 0.02 
and cr = 0.18. The parameters were adjusted so that the underly- 
ing contribution of the RGB stars remains a smooth function of 



the absolute magnitude. RGB stars contribute 49% of all RAVE 
spectra with -1.5 < My < -0.5 and 4500 K< Tes < 5000 K. If 
the absolute magnitude interval is broadened to -3 < Mj < 1 the 
fraction rises to 62%. So their contribution should be accounted 
for when selecting the red clump stars by their effective temper- 
ature, colour or absolute magnitude. 

The 2MASS survey reports also the K magnitudes for the 
same stars. The histogram using the same selection criteria 
(Fig. |6]3) can again be deconvolved into a smooth RGB con- 
tribution and a Gaussian red clump maximum with the peak at 
Mk = -1.61 + 0.02 and cr = 0.19. RGB stars contribute 45% 
of all RAVE spectra with -2.1 < Mj^ < -1.1 and 4500 K< 
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Teff < 5000 K. If the absolute magnitude interval is broadened 
to -3 < Mk < 1 the fraction rises to 60%. 

The position of the red clump peak i s in agreement with 
the va lues frequently used in the literature. I Valentin! & Munaril 
(1201 Ol) used a carefully selected sample of red clump stars 
to estimate Mk = -1.55 an d Mj - -0.87. Similai'ly 
IVan Hels hoecht & GroenewegenI (l2007l) give Mk = -1.57 + 
0.05 and iGroenewegen ( 2008 ) derives M/s: = -1.54 + 0.04. 




Fig. 6. Histogram of absolute J (panel a) and K (panel b) magni- 
tudes for red clump stars. 2MASS J and K apparent magnitudes 
were combined with the distance calculation using the Padova 
isochrone set to derive the absolute J and K magnitudes. Red 
clump stars are assumed to have 4500 K < Tgff < 5000 K. So the 
top panel presents a histogram of Mj values in Fig. |5}l between 
these two temperatures. The histogram representing all stars 
in this temperature interval is decomposed into an assummed 
Gaussian red clump peak and a smooth red giant branch contri- 
bution. For parameter values see text. 



4.3. Formal errors of derived distances 

The errors in our distances arise from four sources. First, the val- 
ues of temperature, gravity and metallicity as derived by RAVE 
have their uncertainties (see Sec. I2.4l i. Moreover, the isochrone 



solution can be uncertain due to intersections of isochrones in 
the parameter space. Also the isochrone positions themselves de- 
pend on assumptions used in their calculation. Finally, we con- 
vert the absolute Mj magnitude into distance assuming a neg- 
ligible interstellar absorption and a reliable measurement of the 
apparent magnitude by the 2MASS survey. The last point has 
been already discussed and generally does not present a major 
concern. We tried to compensate for isochrone differences by 
performing all calculations three times, each time with a differ- 
ent set of isochrones. The uniqueness of the solution has been 
discussed in detail by BIO. This leaves us with studying the in- 
fluence of the errors of measured stellar parameters on the de- 
rived distances. 

Fig.|2]plots the formal errors of the absolute magnitude and 
distance, both as a histogram and as a cumulative plot. The re- 
sults are very similar when using YY and Dartmouth isochrones, 
while the Padova set yields marginally smaller errors. The rea- 
son is a stronger concentration of stars in the red clump region. 
There are virtually no stars with distance errors smaller than 
10%, which is understandable. On the other hand half of the stars 
have their distances determined to better than 33%. These results 
are based on the formal errors obtained from eq.|4]and the errors 
in the parameter described in Sec. 12.41 




(j(d) / d 



Fig. 7. Formal error distribution histogram (left) and cumulative 
plot (right) for absolute J magnitude (top) and distance (bottom). 
Black lines are for YY isochrones, grey ones for Dartmouth, and 
the dashed ones for the Padova ones. The histograms have bin 
widths of 0.05 mag and include 235 064 spectra that have S/N > 
20 and do not show any peculiarities. 

The errors of the stellar parameters vary across the H-R di- 
agram, so it should be the same also with the errors of derived 
distances. The specific structure of the RAVE sample, which has 
a negligible number of very young stars still settling on the main 
sequence and which lacks large numbers of high-mass main- 
sequence stars, implies that the absolute magnitude Mj could 
be used to label the position of the star on the H-R diagram. In 
effect this means that we can use the vertical coordinate in Fig.|5] 
as an independent coordinate for studying distance errors over 
the H-R diagram. Fig.|8]plots the formal scatter of the absolute J 
magnitude for different luminosity classes. The thick line traces 
the median error and the shaded area marks the 20% to 80% 
percentiles of the error distribution. The dashed line shows the 
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Fig. 8. Error in absolute magnitude as a function of absolute J 
magnitude. The thick lines show running median errors for cal- 
culations using YY (a), Dartmouth (b) and Padova isochrones 
(c). The Boxcar for the running median is ±0.5 magnitudes. The 
shaded area shows the 20 to 80 percentile range of the error dis- 
tribution and the dashed line traces the number of objects per 0.1 
mag wide absolute magnitude bin. 



relative number of objects per absolute magnitude bin. It is en- 
couraging to note that the distance errors are lowest for stars on 
the main sequence and in the red clump, i.e. the most frequently 
occurring objects in our sample. On the other hand the errors 
are the largest for stars with Mj ~ 1.0, i.e. for stars that are just 
leaving the main sequence. This is understandable as this is the 
region in which isochrones of different ages intersect in the (tem- 
perature, gravity) plane, making it possible to find a wide variety 
of possible solutions. Another region with large uncertainties is 
that occupied by cool main-sequence stars. Here the errors are 
large on account of uncertain values of the gravity from RAVE 
(see Sec. 12.4b . For some of these stars RAVE measures gravities 
that are too low, so the stars are placed above the main sequence 
(see Fig.|2]i. Thus both the main-sequence and giant solutions en- 
ter eq.|4]and contribute to an increased scatter Fortunately these 
regions of increased distance scatter are sparsely populated in 
our sample. 

Stellar distances are derived using either of the three 
isochrone sets. These sets are not identical, so it is interesting 
to compare the results. Fig.|9]shows the scatter of the three abso- 
lute magnitude values derived for each spectrum. Absolute mag- 
nitudes for the upper main sequence (Teu - 6500 K, log g ^ 4.0) 
and red clump (Tefr - 4750 K, \ogg ^ 2.5) are almost identical 
for all isochrone sets. At the same time these regions have the 
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Fig. 9. Comparison of absolute J magnitudes derived using ei- 
ther YY, Dartmouth, or Padova isochrone sets. Colors label an 
average dispersion in magnitudes for non-peculiar spectra with 
S/N > 20 in a given (temperature, gravity) bin. Solar composi- 
tion isochrones from the YY set with ages between 1 and 1 1 Gyr 
and with a 1 Gyr spacing are overlayed. Metal-poor isochrone 
sets would be shifted towards higher temperatures (see Fig.|2^). 



largest number of stars in our sample (Fig.|5]i. On the other hand 
the most discordant results are for stars in the transition zone be- 
tween the main sequence and giants. Here some isochrone sets 
favour the main-sequence and others the giant solution. The dis- 
persion is somewhat elevated also for stars more luminous than 
the red clump {\ogg ^ 1.9). The origin of this discrepancy lies 
in differences in the physics included in computations of individ- 
ual isochrone sets. Only the Padova set allows for the red clump, 
so the distances computed with this set are to be preferred for 
giants. 

4.4. Errors of derived distances from repeated observations 

Most of the RAVE objects are observed only once. But ~ 10% 
of the observing time is devoted to repeated observations. This 
allows another check of the derived distance errors. Multiple ob- 
servations of the same object are treated independently in our 
calculations, so the results can be directly compared. Note that 
repeated observations allow for a very efficient testing of the dis- 
tance scatter, but they cannot point to possible systematic offsets 
of the mean. We believe that tests using Hipparcos stars with 
known parallaxes and observations of stars in clusters (Sec. |3]l 
give sufficient confidence that any such offsets are insignificant 
for our results. 

The total number of repeated objects in our sample is 19 714 
for which 50 149 spectra have been obtained. This list is limited 
to objects with identical 2MASS identifiers for which RAVE has 
at least two spectra with S/N > 20 and which do not show any 
peculiarities. Fig.fTOlplots the scatter of repeated observations in 
absolute magnitude and distance. The presentation follows that 
of Fig. |7] The results from repeated observations point to a re- 
markable repeatability of our distance estimates. Half of the ob- 
jects show a distance scatter of less than ^11%, and in 80% of 
the cases the distance scatter is smaller than ^31%. Repeated 
observations have a similar distribution in S/N as the other ob- 
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servations, with values which are only slightly higher. So Fig.|2] 
is representative also of the formal errors of repeated observa- 
tions. Comparison of Fig.Qand Fig.[TO]suggests that our formal 
error estimates are rather conservative. 




Fig. 10. Scatter of repeated observations of the same object. 
Scatter distribution histogram (left) and cumulative plot (right) 
for absolute J magnitude (top) and distance (bottom). The num- 
ber of repeated objects with distances calculated with a given set 
of isochrones is 19 714. The presentation follows the scheme of 
Fig.E] 

The scatter of distances derived from repeated observations 
depends on the type of star considered. We plot the scatter of 
repeated observations in Fig. (TT] using the absolute / magnitude 
as its label, in a manner similar to Fig. [8] for the formal errors. 
The morphological properties of the scatter of repeated observa- 
tions is similar to the formal errors, but the values involved are 
much smaller, except for stars with Mj ^ 1 .0, where the scatters 
are similar. The median error for Solar type dwarfs (My ^ 3.0) 
is ^ 0.2 mag. A similarly low scatter is reached also for red 
clump stars when using the Padova isochrones. The dashed lines 
in Fig.fTTImarking the number of repeated objects closely follow 
the distribution of the whole sample (Fig.|8]l. Again they peak in 
the regions of the smallest errors. 

4.5. Distance distribution 

Fig. [12] shows a distribution of distance moduli fi. The results 
for YY and Dartmouth isochrones are very similar. The two 
peaks correspond to dwarfs at distances of ~ 300 pc and gi- 
ants at ~ 1.5 kpc, with some stars being as far away as 10 kpc. 
Calculations using Padova isochrones find more giants and their 
distribution is more peaked. The difference is due to the red- 
clump evolutionary phase which is present only in this isochrone 
set. 

The observed sample, which is limited to the Southern 
Hemisphere, is influenced by the Galactic structure. This causes 
a directional dependence of the derived distances. BIO plot the 
projections of the sample to the Galactic (x,y,z) coordinates (see 
Fig. 10 in BIO). Here we make a more detailed plot (Fig. [13]), 
dividing the sample into 30°-wide slices perpendicular to the 
Galactic plane. The upper set of graphs depicts the general situa- 
tion, and the bottom one is zoomed-in by a factor of 10. The Sun 
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Fig. 11. Scatter of absolute magnitudes calculated for repeated 
observations of the same object as a function of absolute J mag- 
nitude. Thick lines show running median errors for calculations 
using YY (a), Dartmouth (b) and Padova isochrones (c). The 
Boxcar for the running median is ±0.5 magnitudes. The shaded 
area shows the 20 to 80 percentile range of the error distribution 
and the dashed line traces the number of repeated objects per 0. 1 
mag wide absolute magnitude bin. 
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Fig. 12. Distribution of distance moduli (//) for the three 
isochrone sets: YY (black), Dartmouth (grey), and Padova 
(dashed line). The bin size is 0.2 magnitudes. 
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Fig. 13. Galactic vertical cross-sections of the sample for dis- 
tances derived using the Padova isochrone set. Each plot con- 
tains stars with Galactic longitudes within 15° of the labelled 
values. (ipLANE is the cylindrical projected distance along the 
Galactic plane with negative values for stars within the Solar 
circle, z is the distance from the Galactic plane. The density 
contours mark linearly spaced surface densities between 0.001 
and 0.3 stars/pc^. The upper panel shows the whole sample of 
235 064 stars and the bottom panel shows an enlargement of the 
densest parts around the Galactic midplane. 



is located at the origin of the coordinate system, so its distance 
from the Galactic plane is neglected. About half of each graph 
is empty because we are limited to the southern skies. Also the 
general direction of the Southern galactic pole has been observed 
infrequently. This can be seen also in Fig. 21 of Z08. 

The zoomed-in set in Fig. [T3] reveals an expected property. 
Many of the isodensity contours do not follow circular arcs 
around the Solar position. They tend to be horizontal, i.e. par- 
allel to the Galactic plane. This is expected and is a direct con- 
sequence of us living in a disky galaxy. But the fact that we can 
see it in Fig. [T3] supports the notion that the measured distances 
are rather accurate. 

Motivated by this result we plot the histogram of distances 
from the Galactic plane (Fig. [T4l i. again neglecting the Solar 
displacement off the plane. The result seems rather important. 
The RAVE sample generally avoids regions within ~ 20° of the 
Galactic plane. So the histogram values increase until we hit the 
edge of the disk. The decreasing histogram values for stars with 



|z| > 200 pc are well represented by two exponential compo- 
nents, one of them with a very small scale-height. We believe 
such properties illustrate the capacity of RAVE to address the is- 
sues of Galactic structure, using its measurements of kinematics, 
temperature, gravity, metallicity, and now distance. But to draw 
secure conclusions about the properties of the thin and thick 
disks, even if immediately suggestive by the double exponent 
in Fig. [14] requires further study. In particular, a detailed assess- 
ment of possible selection effects and a comparison to models of 
the Galactic structure is needed. Therefore a complete discussion 
is deferred to a separate paper 

4.6. The Catalog 

The final catalog contains distances for 235 064 spectra of 
210.872 different stars. These are RAVE spectra with spectro- 
scopically determined values of radial velocity, surface temper- 
ature, gravity and metallicity. Furthermore these spectra have 
S/N > 20 and were not identified as peculiar by the automated 
pipeline. 

A detailed description of the catalog is given in Table [3] 
Note that distances are given in distance moduli (ji) and not as 
a plain distance id). The reason is a nicer behaviour of error- 
bars. Errors of distance moduli are usually symmetric, while the 
corresponding distance errors are not. The catalog entries will 
be made publicly available through the CDS service alongside 
the regular RAVE public data releases. Currently the latest is 
the second da ta releas e (Z08) with the third one close to pub- 
lication (Sieb ert et al. 1 12010). Distances to RAVE stars are de- 
rived from values of stellar parameters and their adopted errors. 
Should these change in the future we will recalculate the dis- 
tance set and report on the results by updating the corresponding 
astro-ph and CDS entries. 




Fig. 14. Histogram of distances from the Galactic plane if the 
Sun is assumed to be very close to it. The distances are based on 
the Padova isochrone set. 



5. Discussion and conclusions 

The paper uses spectroscopically determined values of stellar pa- 
rameters and the assumption that the star is a single object un- 
dergoing normal stellar evolution to determine its most likely 
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Table 3. Description of the catalog of spectroscopic distances, available at the CDS. 



Column 


Cliaracter 


Format 


Acronym 


Units 


Description 


number 












1 


H7 


17S 


Object-id 




Object name 


2 


19-30 


12.8F 


RA 


degrees 


Right Ascension (J2000.0) 


3 


32-43 


12.8F 


DE 


degrees 


Declination (J2000.0) 


4 


44-50 


6.1F 


RV 


km s"' 


Heliocentric radial velocity 


5 


52-57 


6.1F 


eRV 


kms-' 


Formal RV error 


6 


59-64 


6.1F 


pmRA 


mas yr"' 


Proper motion in RA' 


7 


66-71 


6.1F 


epmRA 


mas yr"' 


Error in pmRA' 


8 


73-78 


6.1F 


pmDE 


mas yr"' 


Proper motion in DE' 


9 


80-85 


6.1F 


epmDE 


mas yr"' 


Error in pmDE' 


10 


87-92 


6.3F 


Jmag 


magnitude 


■/2MASS apparent magnitude 


11 


94-99 


6.3F 


Kmag 


magnitude 


A:2mass apparent magnitude 


12 


101-108 


8S 


Obsdate 


yyyymmdd 


Date of observation 


13 


110-119 


lOS 


Fieldname 




Name of the field 


14 


121-123 


3G 


FibNum 




Fiber number 


15 


125-130 


6G 


Tetf 


K 


Effective temperature 


16 


132-136 


5.2F 


logg 


log(cm s"^) 


Surface gravity 


17 


138-142 


5.2F 


[M/H] 




Metallicity 


18 


144-149 


6.2F 


S/N 




Signal to noise per pixel 


19 


151-157 


7.3F 


IJ{YY) 


magnitude 


Distance modulus (YY isochrones) 


20 


159-165 


7.3F 


e^(YY) 


magnitude 


Formal error in /i(YY) 


21 


167-173 


7.3F 


yu(Dart) 


magnitude 


Distance modulus (Dartmouth isochrones) 


22 


175-181 


7.3F 


e/i(Dart) 


magnitude 


Formal error in p(Dart) 


23 


183-188 


7.3F 


/j(Padova) 


magnitude 


Distance modulus (Padova isochrones) 


24 


190-196 


7.3F 


e/i(Padova) 


magnitude 


Formal error in /i (Padova) 



Notes. If the item is not available its value is set to 9.99E8. 



distance. So the results are expected to be very useful for normal 
stars, while the values for peculiar objects (emission line stars, 
double line spectroscopic binaries) may suffer from systematic 
problems. We note that peculiar objects present a minority in our 
sample, and we also try to detect them using an automated clas- 
sification pipeline. The pipeline is conservative, still it classified 
235 064 out of a total of 251 364 spectra with S/N > 20 as nor- 
mal, i.e. a fraction of 93.5%. Only the distances to such stars are 
published. 

The present effort draws from the extensive experience 
gained during our initial work on the problem (BIO). The 
method is similar but we introduce some refinements. First 
we note that isochrones calculated by different groups are 
not identical. So we give the distances using three recently 
published isochrone sets. The sets are those of Yonsei-Yale 
(Demar queetal. 2004), Da rtmouth (Dotteretal. 2008), and 
Padova (iBertelli et al. Il20()8l) . This leaves it to the discretion of 
the user to choose the isochrone set that is best suited to a par- 
ticular problem. Next we make sure that isochrones are spaced 
linearly in time and locally weighted so that the prevalence of 
low-mass stars is taken into account (we showed that the last 
point is not crucial for the results we get). Also, we take into 
account that the errors of stellar parameter values vary across 
the parameter space. We confirm that the error values published 
in the latest RAVE data release are rather conservative, still we 
prefer to be careful and avoid more optimistic error estimates as 
suggested by the repeated observations. 

The essential differences between our method and that of 
BIO is this. BIO jittered the data by the observational errors and 
sought the model star that most closely matched each jittered 
data point, and then averaged the individual distances obtained 
in this way. We calculate for every model star the probability 
that it could give rise to the data, given the known observational 
errors. Then we use these weights to average the fluxes of the 



model stars and calculate the absolute magnitude of the star. By 
using the closest match, the BIO scheme gives a large weight 
to rapid phases of stellar evolution which sweep the otherwise 
empty parts of the (temperature, gravity) plane. Here we assume 
that stars are observed at a random point of their lives, that less 
massive stars are more numerous, and that higher luminosity 
ones are surveyed in a larger volume of space. 

The results have been tested using two sets of stars for which 
the distances are known a priori. These are Hipparcos stars 
with well determined values of the trigonometric parallax, and 
stars that are secure members of stellar clusters. Distances de- 
termined from RAVE spectra match the astrometric values to 
within i 22%, with a mean difference of less than 3%. This is 
similar to the results already obtained by BIO. Hipparcos stars 
with well determined values of astrometric distance are almost 
exclusively dwarfs. So we need another set to check the dis- 
tances derived for giant stars, which represent about half of the 
RAVE sample. We used members of open and globular clus- 
ters observed by RAVE. We identified spectra of 12 stars that 
are certain members of 7 open clusters and one globular clus- 
ter We note all isochrone sets yield distances which agree with 
cluster distances to within 10% on average. But only the Padova 
isochrone set includes the red-clump phase, so the distances de- 
rived using the Padova set are the preferred ones for giant stars. 
The results for individual stars show a scatter of ~ 20% with re- 
spect to the published cluster distances. This is encouraging, as 
it shows that distances derived spectroscopically agree well with 
other methods up to a range of several kiloparsecs. 

The H-R diagram using the derived distances now appears 
rather clean, with two well defined groups of stars that present 
the majority of the RAVE sample (Fig. |5]l. These are Solar-like 
dwarfs with Teif ~ 6000 K and giants with Tes ~ 4750 K. About 
half of the latter are members of the red clump, with the rest 
being on the red-giant branch. We derive the absolute magnitude 
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of the red clump as Mj = -0.91 ± 0.02 and Mk = -1.61 ± 0.02 
which is consistent with the values from the literature. 

Distance errors have been estimated internally within the 
method and by using repeated observations of the same object. 
The internal estimates appear to be rather conservative and give 
median errors on distance of 33%. The repeatability of derived 
distances is much better: the distances derived from individual 
spectra of the same object match at a level of 11 %. 
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Fig. 15. Comparison of distances (top axis) and distance moduli 
(bottom axis) for a calculation with Padova isochrones if using 
the /iMASS or i^aMASS magnitudes. 



Derived distances are based on apparent magnitudes in the 
■/2MASS band. We assume that interstellar reddening is negligi- 
ble because the vast majority of RAVE stars lie at high Galactic 
latitudes {\b\ > 20°). Nevertheless it is interesting to compare 
the 72MASS results with the ones derived using the /r2MASS band, 
which is even redder and so even less sensitive to reddening. 
So we calculated the distances to the whole sample using the 
^2MASS bandpass and the Padova isochrone set. Fig. [15] com- 
pares the results for both bandpasses. The average distance ratio 
is dx/dj = 0.983 and the scatter is cr{dKldj) = 0.039. Only 5% 
of the spectra have dKJdj < 0.92 and 5% of the spectra have 
dxld] > 1.04. Comparison of our distances and the Hipparcos 
trigonometric parallax determinations (Sec. 13. lb gives similar re- 
sults for both bands, with the scatter in the distance ratio being 
0.25 for the J band and 0.24 for the K band. Also the results for 
stellar clusters are very similar The scatter of spectroscopic dis- 
tances around the cluster distances from the literature is 22% for 
the J band and 23% for the K band. 

One can test whether about 2% smaller distances derived 
from the K instead of the J band are consistent with a decreased 
influence of interstellar reddening at longer wavelengths. From 
relations for extinction for the Solar energy distribution (Aj ^ 
QM7Eb-v and Ak ^ 0.382£b_,/. JFiorucci & Munari I l2003h 
and from Ak - Aj ^ 5 log(0.983) we get Eg-v ^ 0.07. 
This is larger than the values of integrated colour excess to- 
ward t he South Galactic pole which ran ge from Eb-v - 
0.018 (ISchlegel.Finkbeiner&Da"vislll998h to Eb-v = 0.04 
dTee rikorpi 1990). The RAVE sample has esc \b\ - 3.22, where 
b is the Galactic lattitude. So the difference between our average 



color excess and the values at the South Galactic pole can be 
explained with a simple slab model of the Galactic disk. 

We note that a comparison of distances derived using dif- 
ferent passbands can improve our understanding of reddening 
at high Galactic lattitudes. This will be discussed in a separate 
paper, where we shall also address the rare cases of large differ- 
ences in the derived distances which may point to spectral pe- 
culiarities. The photometric effect of interstellar reddening can 
then be confirmed spectroscopically from RAVE data by mea- 
surement of the equivalent width of the diffuse interstellar band 
at 8620 A (iMunari et al. 2008) . We conclude that distances de- 
rived using either the J or the K band are very similar, but a bit 
smaller measurement errors of apparent magnitudes favour the J 
band as our baseline. 

The main result of this paper is a catalog of newly derived 
stellar distances. The paper discusses distances to 235 064 stars. 
The list will increase with time as the RAVE survey progresses. 
The public version of the catalog will be accessible through the 
CDS service, with its size growing along with the public RAVE 
data releases. Note that this is not the only current effort to de- 
rive distances to RAVE stars. Our results are based on the work 
of BIO, and there is also a separate project of distance deriva- 
tion using a Bayesian approach (iBurnett & Binnev 112010 ). We 
believe the availability of complementary computations will be 
of assistance to users, who can then select those best suited to 
their particular problem. 

The present results place RAVE dwarfs at typical dis- 
tances of ~ 300 pc, and giants at 1 or 2 kpc, even 
though there are objects 10 kpc away. In fact RAVE 
reached even the brightest LMC stars (Mun ari et al. I l2009h . 
So the RAVE survey is well placed bet ween the more 
local Geneva Co penhagen Survey (GCS. .Nordstrom et al. I 
12004 iHolmberg. N ordstro m^ Andersen I l2009h and the stel- 
lar component (SEGUE) of the Sloan Digital Sky Survey 
( Abazajia n et al. 2009). Two thirds of the GCS stars are closer 
than 100 pc and 95.5% are closer than 200 pc from the Sun, 
while most of the SEGUE stars are members of the halo. For ex- 
ample the full sample of SDSS/SEGUE DR-7 calibration stars 
contain s only ~ 2400 stars within 1 kpc from the Galactic plane 
(lBeersl l2009). These results make the RAVE survey ideal to 
study the properties of the Galactic thin and thick disks. The 
flatness of the stellar distribution can be seen also directly from 
the shape of the isodensity contours in cross-sections perpen- 
dicular to the galactic plane. This confirms the capability of the 
RAVE survey to address some of the most fundamental ques- 
tions of Galactic astrophysics. The possibility to upgrade the 
RAVE metallicity measurements with individual element abun- 
dances (Boeche et al. 2010) only reinforces such a conclusion. 
Individual applications of this dataset with full 6-dimensional 
kinematical and physical information will be discussed in sepa- 
rate papers. 
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done on the same sample of Hipparcos stars as the one reported 
in Table 1 where the corresponding scatter equals 0.25. A sim- 
ilar test can be performed using members of star clusters. The 
scatter of the spectroscopic to cluster distance ratio is 0.29 for 
the distance calculations using Padova isochrones and reported 
in the main paper The new values of stellar parameters diminish 
the scatter to 0.205. This value is obtained for both choices of er- 
ror estimation: the conservative external error estimates and the 
internal ones obtained from repeated observations. We conclude 
that the best match between our spectroscopic distances and the 
reference values is obtained when using the Padova isochrone 
set, new values of stellar parameters and errors determined from 
repeated observations. 



Note added in proofs 

New values of stellar parameters became available after ac- 
ceptance of the paper. The original results are based on the 
param eter determination pipeline as described in IZ witter et al. I 
(|2008<). while the d erivation of new values is explained in 



ISiebert et al. I (l2010l) . Another recent progress is a redetermi- 
nation of errors of stellar parameters across the H-R dia- 
gram. Instead of th e conservative external eiTors reported in 
IZ witter et al. I ( l2008h we could use an extensive RAVE subset 
where a certain star was observed more than once. Errors on 
stellar parameters were estimated from the scatter of indepen- 
dent parameter determinations for a given star and averag ed over 
stars in the given bin in temperature and gravity ( Siebert et al. I 
I2OIOI) . This internal error seems trustworthy, as there are no ap- 
preciable zero-point offsets between our stellar paramet er values 
and t hose of the reference datasets (see Table 5 in Zwitt er et al. I 
l2008h . The distance moduli have been calculated for the Padova 
isochrone set using the new values of stellar parameters. The 
assumed errors were either the conservative values used in the 
accepted version of the paper or the ones determined from re- 
peated observations. Both distance moduli and their errors add 
four columns to the catalog of spectroscopic distances, available 
via CDS. 

New values of stellar parameters yield distances which are 
generally closer to the reference values than the ones reported 
in the main paper The scatter (cr) of the spectroscopic (d) to 
trigonometric (d^) distance ratio is 0.28 and 0.22 if the external 
and internal error estimates are used, respectively. This test was 
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